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EFFECTS OF SOLAR ACTIVE REGIONS ON MERIDIONAL FLOWS 
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ABSTRACT 

The aim of this paper is to extend our previous study of the solar-cycle variations of the meridional flows 
and to investigate their latitudinal and longitudinal structure in the subphotospheric layer, especially their vari- 
ations in magnetic regions. Helioseismology observations indicate that mass flows around active regions are 
dominated by inflows into those regions. On average, those local flows are more important around leading 
magnetic polarities of active regions than around the following polarities, and depend on the evolutionary stage 
of particular active regions. We present a statistical study based on MDI/SOHO observations of 1996-2002 
and show that this effect explains a significant part of the cyclic change of meridional flows in near-equatorial 
regions, but not at higher latitudes. A different mechanism driving solar-cycle variations of the meridional flow 
probably operates. 

Subject headings: Sun: atmospheric motions - Sun: magnetic fields - Sun: activity 



1. INTRODUCTION 

Meridional flows are one of the integral properties of the 
solar convective envelope dynamics. They are essential for 
the magnetic flux transport from the activity belt towards the 
solar poles and ther efore play a signi ficant role in the polar 
field reversals (e.g. IWang et al.l [19891) . The average merid- 
ional flow profile (MFP) derived by local helioseismology 
is often used as an observational input in models describ- 
ing the solar dynamo an d magnetic flux transportation (e.g. 
iDikpati & Gilmanll2006h . As we pointed out in our previous 
study (ISvanda et al.ll2007h . it seemed that the flows around 
large active regions affect the mean MFP obtained as a lon- 
gitudinal average of the North-South component of the flow 
velocity map. In this paper we investigate this effect in more 
detail and show that it may explain most of the observed vari- 
ations of the meridional flow speed in low-latitude region, but 
not at high latitudes. 

2. data AND analysis 

Local helioseismology methods provide an unique in- 
formation about the subsurface dynamics of the Sun (e.g. 
iHaber et^l. 200 41). In par ticular, the time-distance helioseis- 
mology (Duvall et al.lll993 ) is a tool that allows to use inver- 
sions of solar oscillation observat ions to infer the s tructure 
and topology of subsurface flows (lKosovichevlll9 96^. Solar 
acoustic waves (p modes) are believed to be excited in the 
upper convection zone. They travel between various points 
on the surface through the interior and are perturbed by mass 
flows along the path of propagation. The mass flow velocities 
in the interior are inferred from the differences of travel times 
from the central point to the surrounding annuli and the travel 
times from the surrounding annuli to the central point. 

The travel-time inversions are applied to full-resolution 
full-disk 1-min -cadence data observ ed by Michelson Doppler 
Imager (MDI; IScherrer et all 1 19951) on-board the Solai" and 
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Heliospheric Observatory (SoHO). These data are available 
for only approximately two months a year uninterruptedly. 
From those data seven non-consecutive synoptic flow maps 
covering seven solar rotations during the period of 1996-2002 
were constructed (Zhao & Kosovichev 2004). The data from 
the declining phase of the solar cycle were not used because 
they contain some probably instrumental issues and we de- 
cided not to use them for this particular study. For study of 
the relationship between the flows and magnetic fields synop- 
tic maps of the line-of-sight component of magnetic field were 
constructed using the same algorithm based on MDI magne- 
tograms. 

In seven processed Carrington rotations we selected 92 
magnetic active regions to investigate the influence of the 
flows in and around them to the total longitudinally averaged 
meridional flows by applying masks to isolated local flow pat- 
terns. For a proper masking, we need to determine the area of 
influence of active regions. We started with a rectangular box 
around each selected region, in which the magnitude of the 
line-of-sight component of magnetic field is above 25 Gauss. 
In this box we calculated a characteristic perturbation 5v, 
defined as a mean of the absolute differences of meridional 
speed with respect to MFP calculated separately of each Car- 
rington rotation, when the box around particular active region 
was masked. Then we expanded the box around this region 
while calculating the deviations 5v in this box. The expanded 
box, where 5v reaches its first maximum, we considered as 
an area of influence of particular active region. Other max- 
ima are attributed to the presence of other magnetic region, 
which appeared in the expanded box. The regions of influ- 
ence around each magnetic area were masked in the synoptic 
map. The size of the region of influence was on average twice 
as large as the initial 25-Gauss box. 

3. RESULTS 

Variations caused by active regions — The new analysis con- 
firmed the results of our previous study (ISvanda et al.ll2007h . 
The local flows around strong active regions affect the mean 
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Fig. 1 . — The longitudinally averaged meridional component of the time- 
distance helioseismology flow map at 3^.5 Mm depth (solid line) and the 
same quantity when active regions are masked (dotted line) for Carrington 
rotation No. 1975 (April 2001). Error-bars have size of ~2 ms"'. 

meridional flow obtained as a longitudinal average of the 
North-South component of the synoptic flow map. On av- 
erage, the large-scale inflow pattern in the activit y belt pre- 
vails, which is in agreement with fi ndin gs by Basu & Antial 
(I2003h. IZhao & KosovichevI (|2004 . iHaber et al.. (20041), and 
iKomm et al.l (l2006h . For example, in Fig. [T] one can clearly 
see that MFP calculated excluding active regions (henceforth 
called as the reference meridional flow profile - RMFP) varies 
less with latitude than when the flows associated with active 
regions are included. The deviations between the profiles with 
and without the active region flows are up to 5 m s"' . On av- 
erage, the positive deviations in the areas of influence prevail 
in the Northern hemisphere (by 0.08 ms"'), while the nega- 
tive ones prevail in the Southern hemisphere (by 0.10 ms~'). 
This may be interpreted that the inflows into active regions 
are stronger on their equatorial side than on poleward one, in 
both hemispheres. Such small differences are not statistically 
significant; thus to confirm this, a detailed stud y is needed . 
An analytical model of torsional oscillations by S pruit! (l2003l) 
predicted such inflows into the activity belt with characteristic 
amplitude of 6 ms~' . Simila r value s were predicted also by 
numerical model by iRempell (l2006t) . Both models assumed 
a mechanism of thermal forcing for the low-latitude torsional 
oscillations branches. The inflow in the activity belt was a 
side-effect of these models, but it resembles the observed in- 
flow in active regions. 

We investigated a possible dependence of the characteristic 
perturbation of the flows at different depths in individual ac- 
tive regions on the amount of magnetic flux present in their 
areas. As a measure of the characteristic perturbation Jvmer 
we used a mean of the absolute differences between the mea- 
sured meridional component of the flow field and RMFP at 
given latitude calculated over the region of interest. 

The dependence of Jv^er on the total unsigned flux for all 
92 active regions is shown in Fig. |2] Other characteristics, 
such as a geometrical average of deviations over the whole 
active region, display a similar dependence. The perturba- 
tions become less significant with the increase of depth in the 
convection zone. Fig.|2]indicates that the active regions con- 
taining more magnetic flux cause more significant deviations 
from RMFP than those containing less flux. 

Examples of synoptic maps of the meridional flow com- 
ponent for Carrington rotations CR 1911 (July 1996) and 
CR 1975 (April 2001) are shown in Fig. [3] These maps con- 
tain a wide range of local motions of various scales, but their 
structure is far from the idealized single-cell picture. In the 
map of CR 1975 one can see very large and complex active 
regions with converging flows (inflows) prevailing in the lead- 



ing parts of the active regions, while the trailing parts do not 
seem to differ in topology of flows from the quiet Sun regime. 
This e ffect was already suggested by Kosovichev & Duvalll 
(|2006|). In those examples we see that converging flows ex- 
tend from the active region far in the poleward direction and 
form a surge-shaped structure. 

Asymmetry between the flows around leading and following polari- 
ties of active regions — The question is whether this asymmetry 
between the flow around the leading and following polarities 
is common for all active regions. To study this issue we di- 
vide the box around each active region in two parts separating 
the regions of the leading and following polarities. As a mea- 
sure of the flow deviations in leading and trailing polarities 
we use the difference A between Sv^^^ in both polarities of 
each active region, i.e. 

A = (5viner(leading)-(5vn,er(ti-ailing) . (1) 

We do not find a clear correlation of A with any characteristic 
of magnetic field of the active regions. However, a histogram 
of A (Fig. nil shows that the leading and following polarities 
differ quite often as the positive part contains more cases than 
the negative one. 

The distribution of A seems to be bimodal. The second 
peak (around 0.4 ms"') suggests an existence of group of ac- 
tive regions, where inflows in their leading polarities are sig- 
nificantly stronger than in trailing ones. Some examples of ac- 
tive regions, for which the flows have the behavior described 
above, can be seen in Fig. [3] There is no typical characteristic 
or feature among the active regions responsible for the sec- 
ond peak in the distribution of A in our sample, except for the 
tilt angle, calculated from positions of the centres of mass of 
both polarities. The tilt angle is on average larger (38 degrees 
vs. 19 degrees) for the active regions that form the second 
peak in the distribution of A. However, those numbers are 
not convincing due to the small sample. Perhaps, this effect 
ca n be related to th e age of active regions. As pointed out, e.g. 
by iHowardI (1 1 9941) . the tilt angle of bipolar sunspot groups on 
average may increase during their evolution. 

Our visual inspection shows that most of 1 1 active regions 
making the second peak of the histogram have a trailing po- 
larity dispersed probably by diffusion, while the leading one 
seems to be more compact. A quick look on SolarMonitor.org 
and Mt. Wilson drawings archive showed that these active re- 
gions are old and highly evolved with complex topology. A 
detailed investigation of this effect requires a larger dataset of 
helioseismic and magnetic observations of active regions. 

Latitudinal variations of the meridional flows — In Fig. |5] the 
evolution of the RMFP averaged over depths of 0-9 Mm with 
solar-cycle is given, from which the profile in minimum of 
solar activity (CR191 1) was subtracted. We basically see that 
the meridional flow slows down at higher latitudes with the in- 
creasing magnetic activity, and that variations are larger when 
the magnetic regions are included in the synoptic flow maps. 
The differences between corresponding curves in Figs.|5^ and 
I5J5 are statistically significant. In Fig.|5}; the variations caused 
by the magnetic regions averaged over two latitudinal belts in 
both hemispheres are given. The opposite sign of the merid- 
ional flow in both hemisphere was taken into account. We see 
that the deviations from RMFP caused by the active regions 
increase with activity cycle, which is consistent with results 
described in previous sections. 

A large amount of cycle-related variation can be ex- 
plained by flows around the active regions in near-equatorial 
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Fig. 2. — A dependence of the characteristic variations of the meridional flow on the total unsigned flux localized in active regions for the flows at the depth of 
0-3 Mm (left) and 3^.5 Mm (right). The dependence in both depth intervals is very similar, but the magnitude of deviations is smaller at greater depth. The '+' 
signs are for areas in the Northern hemisphere, ' X ' sings for the Southern hemisphere ones. 
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Fig. 3. — Examples of synoptic maps of the meridional (North-South) component at the depth of 3^.5 Mm smoothed and binned by 10 heliospheric degrees. 
The positive meridional speed is Northward. The magnetic active regions are displayed with contours: white contours show the positive polarity, the black ones 
show the negative polarity. These maps demonstrate a longitudinal structure of meridional flows close to minimum (CR 1911) and maximum (CR 1975) of solar 
activity. The en'ors of individual flow velocity measurements are ~2 ms"' aZhao & Kosovichev|[20Q4l) . 

belt. If we quantify this effect calculating value i9 = (vmer - 

Vmer,ref)/ (Vma- Vmer,CRI91 1 ), whcrC V^er is MFP including floWS 

around the active regions, Vmer.ref the one with active regions 
masked, and Vmer.CRigii is MFP in minimum of solar activ- 
ity, we may estimate the importance of the flows around the 
active regions on the cycle-related variations of MFP. For 
the equatorial region (latitudes of to 15 degrees) we find 
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Fig. 4. — Histogram of the characteristic deviations of meridional compo- 
nents from the RMFP. showing the difference between the leading and trailing 
polarities of active regions in our sample. The histogram reflected around AO 
is show in dotted line to enhance the visibility of the asymmetry. 



d = (0.65 ±0.21) in the Southern and = (0.34 ±0.26) in 
the Northern hemisphere, and for higher latitudes (15 to 40 
degrees) we find = (0.20 ± 0.15) for the Southern and -d = 
(0.37 ± 0.29) for the Northern hemisphere. The North-South 
asymmetry is statistically significant, the effect is stronger in 
the Southern hemisphere, where also the surface magnetic ac- 
tivity is stronger than on the Northern hemisphere. On av- 
erage, the effect of the flows around the active regions is 
not sufficient to explain cycle-related changes of MFP, es- 
pecially at higher latitudes, where remaining cycle-related 
variations may be explained by formation of persistent high- 
latitude counter-cells providing the mass inflow on poleward 
side of activity belt. The peak speed in such counter-cell could 
be estimated by ^10 ms"' in the Northern hemisphere and 
~25 ms"' in the Southern hemisphere. Perhaps, this may be 
due to the existence of deep magnetic field, which could also 
depict North-South asymmetry. We think that these results 
represent an interesting challenge for theories and numerical 
simulations of the solar meridional circulation. 

Longitudinal variations of RMFP — The question whether or 
not the MFP without active regions is the same at different 
longitudes naturally arises. Unfortunately, we cannot provide 
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Fig. 5. — a) - Solar-cycle variations of MFP averaged over depths of 0- 
9 Mm after the subtraction of MFP in CR1911. b) - The same with active 
regions masked in the maps, c) - The difference of a) and b) averaged over 
two latitudinal belts to demonstrate the effect of the active regions on longitu- 
dinally averaged MFP at various phases of the magnetic activity (the negative 
values resemble equatorward flow). 

a clear answer based on the available data. As we can see 
in Fig. |3] the synoptic maps of the large-scale North-South 
flows component show a lot of variations, which are not per- 
sistent for the whole synoptic map. These local disturbances 
of the meridional component of the flow have greater magni- 
tude than the mean flow. The magnitude of the disturbances 
is decreasing with depth. One needs to keep in mind that in 
general, it is impossible to distinguish between the temporal 
evolution of the meridional flow and its longitudinal structure 
in synoptic map, which is composed of daily strips around 
the central meridian. However, the assumption that persistent 
component of meridional flow evolves slowly with character- 
istic time of several days seems reasonable, which justifies 
our ap proach . It wa s described in numerical simulations (e.g. 
iMiesch et al.ll2008l) that meridional flow has a strong fluctuat- 
ing component, magnitude of which is larger than magnitude 
of the mean temporally averaged meridional flow. Some fluc- 
tuations surviving for many hours are also detected in our data 



and the mean poleward meridional flow appears only after av- 
eraging. Perhaps, similar fluctuations happen in the real Sun 
as well. 



4. CONCLUSIONS 

From our study it becomes clear that the large-scale mass 
flows around active regions significantly influence the mean 
meridional flow. This is an important conclusion for the 
dynamo and solar cycle models, which use the longitudi- 
nally averaged meridional flow from local helioseismology 
as a measure of the mean meridional flow. Profiles calcu- 
lated simply by taking a longitudinal average do not repre- 
sent the mean meridional circulation in and between activity 
belts. For particular active region, the deviations caused by 
the circulation flows around them depend on the total mag- 
netic flux in the area. On average, the local circulation flows 
are represented by inflows in the activity belt, as predicted by 
Spr uit (2003) and obseryed by i.e. Basu & Antia (2003) or 
Zhao & Kosovichevl(f20l34l) . 

Moreover, it seems that, on average, the local flows are 
more important in the areas of leading magnetic polarities 
of active regions. There is a group of active regions in our 
sample, where local flows are significantly more important in 
the leading polarities. This effect may be related to the evo- 
lution of active regions and requires a detailed study based 
on extended dataset describing three-dimensional structure of 
flows in and around active regions in various stages of their 
evolution. Those deviations explain a significant part of the 
cyclic change of the meridional profile shape in the equato- 
rial region, but not at higher latitudes. A different mecha- 
nism must operate to explain all cyclic variations. Using the 
currently available dataset we find an evidence that the latitu- 
dinal variations of the meridional flow at high latitudes may 
be related to subsurface magnetic activity, but at present state 
this is inconclusive. Another interesting problem is how the 
asymmetry in the strength of the local inflows around leading 
and following polarities of active regions affects the magnetic 
flux transport to the polar regions and Sun's polarity reversals. 
As we pointed our in Svanda et al.l (l2007i) . this phenonenon is 
not important for the meridional flux transport averaged over 
a few Carrington rotations. We expect that in detail this con- 
clusion may be changed. This problem can be studied when 
more local helioseismology data will become available from 
Solar Dynamics Observatory. 
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